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Abstract 

We present a Lattice QCD determination of the chiral quark condensate based 
on a new method. We extract the quark condensate from the operator product 
expansion of the quark propagator at short euclidean distances, where it repre¬ 
sents the leading contribution in the chiral limit. From this study we obtain 
{qq}^^(2 GeV) = —(265 ± 5 ± 22 MeV)^, in good agreement with determinations of 
this quantity based on different approaches. The simulation is performed by using 
the 0(a)-improved Wilson action at /3 = 6.45 on a volume 32^ x 70 in the quenched 
approximation. 


1 Introduction 


An accurate determination of the chiral quark condensate is a task of prime interest. Its 
non-vanishing value signals the spontaneous breaking of chiral symmetry in QCD and, 
quantitatively, it is related to the pseudo-Goldstone bosons mass spectrum. 

Due to the purely non-perturbative nature of the quark condensate, its estimate is 
rather challenging. Traditional approaches have been based on QCD sum rules (a review 
of these techniques can be found in refs. IB 121). In the last years, hrst principle deter¬ 
minations of the quark condensate have been provided by Lattice QCD calculations, and 
the accuracy of these results is expected to systematically improve in time. The standard 
method to extract the quark condensate from lattice calculations exploits the well known 
CMOR formula [21-[HI- Alternative techniques have been also investigated, based on the 
e-expansion of QCD in a small volume 0-i and on the study of the Coldstone pole 
contribution to the pseudoscalar quark Creen function 

In this paper, we present an exploratory Lattice QCD determination of the chiral 
quark condensate based on a new method. We study the quark propagator in coordinate 
space and its operator product expansion (OPE) |T2] at short euclidean distances. The 
OPE is a powerful technique that systematically includes non-perturbative corrections 
and parameterizes the non-trivial properties of the QCD vacuum in terms of conden¬ 
sates uni. We extract the quark condensate {'ijj'ip) by evaluating the quark propagator at 
short distances on the lattice, and comparing the result with the OPE prediction, 

~ + ... ( 1 ) 

(tZ' ) tZ' 

where the dots represent higher powers of x"^ and of the quark mass m} Our hnal result 
for the chiral quark condensate, renormalized in the MS scheme at the scale /i = 2 GeV, 
is _ 

(gg)^®(2 GeV) = -(265 ± 5 ± 22 MeV)^ , (2) 

where the hrst error is statistical and the second systematic. This result is in good 
agreement with those obtained from the other methods listed above. It also provides a 
remarkable non-perturbative test of the OPE predictions at short distance in QCD. 

The OPE of the quark propagator can be also performed in momentum space, from 
which a determination of the quark condensate might be possible as well. When working 
on the lattice with Wilson fermions, however, the leading contribution to the OPE in 
momentum space is a constant term induced by discretization effects. Though vanishing 
in the continuum limit, this term is dominant at hxed lattice spacing with respect to the 
mass and the condensate contributions, whose coefficients are suppressed by 1/p^ and 1/p^ 
respectively [Hj. In coordinate space this major obstacle is bypassed, since the Fourier 
transform of the unphysical term is a discretized delta function, whose effect is negligible 
at distances larger than few lattice spacings. 

Another advantage of the approach studied in this paper is that it greatly simplihes the 
renormalization procedure. Specihcally, once the quark propagator on the l.h.s. of eq. m 

^Throughout this paper we use the notation x = 


1 




is renormalized, all contributions appearing on the r.h.s. turn out to be expressed in terms 
of renormalized quantities. In particular, the determination of the chiral quark condensate 
in this approach does not require the evaluation of the corresponding renormalization 
constant. 

The applicability of the OPE to correlation functions evaluated on the lattice at hxed 
value of the lattice spacing a relies on the existence of a short distance region where the 
conditions 

X ^ 1/Aqcd (3) 

are both satished. The upper bound in eq. guarantees that the Wilson coefficients 
entering the OPE at the typical scale fj, = 1/x can be evaluated in perturbation theory. 
The lower bound must be satished in order to keep under control discretization effects. 
In the present study, though we use an (P(a)-improved action and the value of the inverse 
lattice spacing is as large as a~^ ~ 4 GeV, we hnd that in the region x ^ 1/Aqcd 
discretization effects in the quark propagator are not negligible. These effects are in fact 
responsible for most of the systematic uncertainty quoted in eq. In order to reduce 
their contribution, we have followed a procedure similar to the one applied in ref. ng: 
we have corrected the lattice results for the quark propagator by the lattice artifacts 
computed in the free theory, thus reducing their size from 0{a^) to O(asa^). 

We now summarize the procedure followed in this study and present the plan of the 
paper. 

- In sect.2, we derive the OPE of the quark propagator in coordinate space, by including 
QCD corrections up to the next-to-leading order (NLO). 

- Details of the lattice simulation are presented in sect.3, where the tree-level correction 
of lattice artifacts is also discussed. 

- In sect.4 we compute the renormalization constant of the quark propagator non-pertur- 
batively in the X-space scheme. The X-space method has been proposed in ref. cni, and 
applied in to compute the renormalization constants of bilinear quark operators. Our 
result for the quark held renormalization constant, converted to the MS scheme, reads 

= 2 GeV) = 0.871 ± 0.003 ± 0.020 , (4) 

in good agreement with the result obtained in ref. Hg by using the non-perturbative 
RI-MOM method. 

- In sect.5 we evaluate the chiral quark condensate by htting in coordinate space the quark 
propagator, extrapolated to the chiral limit, to its OPE. A second estimate is obtained by 
hrst using the OPE at hnite values of the quark mass and then extrapolating the result to 
the chiral limit. Different functional forms are considered in the hts, and the differences 
among the results are taken into account in the estimate of the systematic error. The two 
approaches give completely consistent results. 

- The hnal result quoted in eq. m is presented in sect.6, where we discuss in details the 
evaluation of the systematic error. 

- Finally, we sketch in the appendix the NLO QGD calculation of the Wilson coefficients 
entering in eq. (|H) 
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2 OPE of the quark propagator in coordinate space 


The quark propagator can be expressed in terms of two scalar form factors, Si(a;) and 
T, 2 {x), which are defined from 



The leading terms in the OPE of Ei(a;) and T, 2 {x) can be read from eq. O: 

^ 2 ix) = Cm{x) m - ^ C^^{x) x"^ + ■ ■ ■ (6) 

where, at variance with eq. (HD, the Wilson coefficients Cj{x), Cm(x) and C^^{x) are 
normalized to unity in the free theory. is the number of colors and the quark condensate 
is defined as 

(#) = (C(O)C(O)). (7) 

where a summation over repeated color and spin indices is understood. 

By using the known two-loop results for the quark field and the quark mass anomalous 
dimensions in QCD, a simple one-loop calculation provides the renormalization group 
improved expressions for the Wilson coefficients in eq. ( 0 , at the NLO. The main steps 
of the calculation are given in the appendix. We find, in the MS scheme. 


^MS 




2'k‘^ 


l/x) 


l-2^i^CK(7E-log2) 


( 8 ) 


^MS 


2 1/x) 

f 1 


l-2^^^C^e(7E-log2) 


1 

"4W 


Wn(/U, l/x) 


1 + 


dvr 
as{l/x) 




1 + 2 


dTT 
as{l/x) 
dTT 


C',7(d-2(e-3)(7i,-log2)) 




(9) 


C'f (1 — ('C + 3)(ye — log2)) 




The terms in square brackets represent the Wilson coefficients at the scale fi = l/x, 
whereas l/x), with I = 'ip,m, are the NLO evolution functions. 


Wi{n, l/x) 
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( 10 ) 


The coefficients of the beta function and of the quark mass and quark field anomalous 
dimensions at the LO and NLO read: 


Po — 


11W - 2n 


f 


Pi — - 2CFnf 


3 



























K 

0.1349 0.1351 0.1352 0.1353 

am^,(2 GeV) 
GeV) 

0.0305(4) 0.0227(3) 0.0188(2) 0.0149(2) 

0.0288(3) 0.0215(2) 0.0178(2) 0.0141(1) 


Table 1: Quark masses in lattice units for the non-perturbatively O{a)-improved Wilson 
action at f3 = 6.45. The results are taken from ref. (where they are quoted in the 
RI-MOM scheme at the scale /i = 3 GeV ). 

/Q7 in \ 

7m = SCf , 7^ = Gf iVc + - —77./j (11) 

7? = -2fCf. , 7‘ + 

where Cp = (iV^ — l)/{2Nc), ^ is the gauge parameter (^ = 0 in the Landau gauge) and 
uj is the number of active flavors (n/ = 0 in the quenched approximation). 

The result in coordinate space for the Wilson coefficient of the quark condensate given 
in eq. m corresponds to the one obtained in ref. ffi momentum space. Eqs. (jHl) and 
(jni) will be used in sects.4 and 5 to extract the quark held renormalization constant and 
the chiral quark condensate with NLO accuracy in the MS scheme. 

3 Analysis of discretization effects 

In this section we present the details of the lattice simulation, illustrate the results ob¬ 
tained for the bare quark propagator and discuss the free theory correction implemented 
in order to reduce the lattice artifacts. 

We have generated 180 gauge conhgurations in the quenched approximation with the 
non-perturbatively (!l(a)-improved Wilson action on a volume 32^ x 70 at /3 = 6.45. As a 
value of the inverse lattice spacing we use a~^ = 3.87(19) GeV, as obtained in ref. [1^] from 
the studies of the quark-antiquark potential m and by using in input the reference scale 
a~^{j3 = 6.0) = 2.0(1) GeV.^ We have computed the quark propagator at four values of 
the hopping parameter, k = 0.1349, 0.1351, 0.1352, 0.1353, corresponding to light quark 
masses in the range mg/2 ^ m ^ m^. The corresponding values of the renormalized quark 
masses have been obtained in ref. Hi from the study of both the vector and axial-vector 
Ward identities, and are given in lattice units, in the MS scheme, in table [H These values 
have been used in the study of the OPE of the quark propagator and to perform the chiral 
extrapolations of the quantities we are interested in. The statistical errors quoted in this 
paper have been evaluated with the jackknife technique. 

We have hxed the Landau gauge on the lattice by minimizing the quantity: 

9 = [Af,Af,{x) A^Ai,{x)] ( 12 ) 

* X 

^Had we used in input rg = 0.5 fm we would have obtained a~^ ~ 4.10 GeV. 
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where V is the lattice volume and is the discretized version of the gauge held 

divergence d^A^. We have required 6 < 5.0 x 10“^ for all the conhgurations used in this 
study. 

In order to compute the quark condensate and the quark held renormalization con¬ 
stant, we have extracted from the quark propagator in the Landau gauge the bare form 
factors Si(x) and S 2 (a;) dehned in eq. ©• The results are shown in hg. mtop) for 
K = 0.1349 as functions of (x/a)^. 

Lorentz invariance requires that, when approaching the continuum limit, the form 
factors should become functions of only. At hxed value of the lattice spacing, however, 
the plots in hg. IH^top) show that points corresponding to the same value of (x/a)^ are 
signihcantly spreaded out. This is true especially at short distances ((x/a)^ ^ 10), where 
discretization ehects are expected to be larger. A better understanding of these ehects 
can be obtained by studying the lattice quark propagator in the free theory. Indeed, in 
the short distance region which is relevant for the X-space method the interacting theory 
is expected to approach the asymptotic free regime, up to small perturbative corrections. 
One hnds that Ei(x) and S 2 (x), computed on the lattice in the free theory, present similar 
deviations from the expected continuum behavior. The free theory results, obtained at 
inhnite volume and in the chiral limit, are shown in hg. ^center), and they can be 
compared with the lattice results shown in the top panels. This similarity suggests that 
one can identify the discretization patterns in the free case in order to subtract them in 
the interacting case of interest ini. The practical implementation of this approach passes 
through the dehnition of the “corrected” form factors. 

In the case of Si(x) we dehne 

( ycont \ 

(13) 

where S5°j*ee(^) ^^^1 S^“j^gg(x) are the free theory form factors computed respectively in 
the continuum and on the lattice, at inhnite volume and in the chiral limit. For hnite 
values of the lattice spacing, the diherence of the ratio 

is a measure of tree-level discretization errors. After the correction of eq. m, we expect 
these errors to be reduced from O(a^) to O(asa^). 

Concerning S 2 (x) one observes that, in the continuum and in the chiral limit, the 
form factor vanishes at any order of perturbation theory. Therefore, in the case of S 2 we 
implement the following correction: 

E—(x) = S2(x)-SJVgg(x), (14) 

where represents a pure discretization ehect. After eqs. (jl3|l and (ED have been 

implemented, we also average the results for the form factors Ei(x) and E 2 (x) obtained 
at lattice points which correspond to the same value of x^. 

The remarkable ehect of the correction on the two form factors is shown in hg. Q 
(bottom). In the following analysis, otherwise indicated, we will always use the corrected 
form factors dehned in eqs. ED and ED- 
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Figure 1: The bare form factors Si (left panels) and S 2 (right panels). We show, from 
top to bottom: the form factors in the interacting theory, as obtained from the lattice 
simulation at k=0.1349; the form factors computed in the free lattice theory, at infinite 
volume and in the chiral limit; the “corrected” form factors defined in egs. m and a- 
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4 Renormalization of the quark propagator in the X- 
scheme 


In this section, we define the X-space renormalization scheme HSI for the qnark propagator 
and discuss the determination of the corresponding renormalization constant. 

The quark field renormalization constant in the Landau gauge X-scheme, is 

determined non-perturbatively by imposing the condition 


= l/x)^i{x) 


5=0 
m—^■O 


s^cont 

^IJree 


(x) 


(15) 


where the value of the form factor in the free continuum theory and in the chiral limit 
is E“j*ee(a;) = l/(27r^). The limit m ^ 0 in eq. guarantees a mass-independent 
definition of the renormalization scheme. It also guarantees that, when using the 0{a)- 
improved Wilson action, the renormalization constant computed from eq. m is auto¬ 
matically (9(a)-improved, without need of further improving the quark field [T7j.^ 

In order to extrapolate the form factor Ei(a;) to the chiral limit we have assumed a 
linear dependence on the quark mass. This dependence describes well the lattice data as 
can be seen from fig. | 2 l where the linear fit is shown for three values of in the range 
of interest. A quadratic fit has been also performed in order to evaluate the systematic 
error involved in the chiral extrapolation. 

By combining the renormalization condition (USD with the NLO evolution function of 
given in eq. (Uni) and considering that the LO anomalous dimension of the quark field 
vanishes in the Landau gauge, one finds at the NLO 


1/a:) Z^{l/x) 


/ _ Yt a.W-a.(l/x) \ f Si%Jxn 

V 2/3(1 4/r A ] 


(16) 


We also note that, in the Landau gauge, the equality of the renormalized form factor Si(x) 
at one loop in the MS and X schemes implies that the NLO anomalous dimensions 7 ^ are 
also equal in the two schemes. In the numerical analysis, as{l/x) has been evaluated at the 
NLO in the MS scheme by using n/ = 0 and the quenched estimate = 0.225(21) GeV 
(obtained from roAg-^L^ = 0.602(48) j2I] and using tq = 0.525(25) fm). 

As already discussed in the introduction, the X-space non-perturbative renormaliza¬ 
tion approach relies on the existence of a window a ^ x ^ 1/Aqcd which permits matching 
the lattice results with the perturbative ones and, at the same time, to avoid the region at 
very short distances affected by contact terms and large discretization effects. In practice, 
in the present study, we consider this condition satisfied in the range 9 ^ (x/a)^ ^ 25 
(the upper bound corresponds to x~^ ~ 1 GeV). 

^In the definition of the C>(a)-improved quark field, 

qi{x) = q{x) + a bqmq{x) + a c'q{Ip + m) q{x) + a CNGi^q{x) 

the second term vanishes in the chiral limit, while the third one produces in the quark propagator a 
contact term in x = 0. The contribution to the quark propagator of the last non gauge-invariant term 
has been found to be practically indistinguishable from the contact term proportional to c ' m 
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Figure 2: Quark mass dependence of T,i{x) at three different values of x^. The dashed 
lines represent the result of a linear fit. 



(x/a)^-range 

fit 

0.871 ±0.003 

[9,25] 

constant 

0.852 ±0.003 

[9,40] 

linear 


Table 2: in the X-space scheme at p. = 2 GeV as obtained from either a constant 

or a linear fit in x^. The quoted errors are statistical only. 

The results for Z^{p = 2 GeV) as obtained from eq. (fTT)|l at different values of x"^ are 
shown in fig. El One observes that even in the region [x/afi = [9, 25] the data show some 
spread, at the level of few percent. This spread is due to discretization errors which remain 
after the free theory correction has been implemented. It represents the main source of 
systematic uncertainty in the evaluation of A second source of uncertainty is due 
the fact that one cannot exclude, even in the fitting region {x/afi = [9,25], a systematic 
dependence of the data on which could be due to higher order contributions to the 
OPE of Ei(a;) neglected in eq. (P). In order to evaluate this systematics, we have evaluated 
Z^(/r) from both a constant and linear fit in {x/afi, by considering in the latter case larger 
intervals in x^ (np to {x/afi = 40). The results of the fits are presented in tabled 

The other sonrces of systematic effects, as those deriving from the determination of 
the lattice scale, the estimate of Aqcd, the difference between linear and qnadratic chiral 
extrapolations and the use of the vector or the axial-vector definitions of the qnark masses 
in these extrapolations, are fonnd to be negligible. As final estimate of Z^{p) we thus 
quote 

Zj(2 GeV) = 0.871 ±0.003 ±0.020 (17) 

where the first error is statistical and the second systematic. The central valne in eq. (ED 
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Figure 3: Values of = 2 GeV) as obtained from eq. HE) for different values of 

(x/a)'^. The solid lines indicate the results obtained from a constant fit in and the 
estimated systematic error (eq. The dashed vertical lines show the range in {x/afi 

where the constant fit is performed. The result of a linear fit in is also shown (dashed 
line). In the latter case, the estimate of = 2 GeV) is given by the intercept. 

is the one obtained from the constant fit in the shorter distance range {x/afi = [9, 25], 
where the contribution of higher power corrections is more suppressed. The few percent 
error on the value of introduces an uncertainty in the estimate of the chiral quark 
condensate discussed in the next section which is completely negligible. 

The vanishing of the one-loop contribution to the form factor Ei(x) in the Landau 
gauge implies that the quark held renormalization constant at the NLO is equal in several 
commonly used renormalization schemes. In particular, 

= Zf{tv) = (18) 

at the NLO. The result in eq. (El can be therefore directly compared to the value 
^ri-mom (2 QgY) = 0.865 ± 0.003 obtained non-perturbatively in ref. [T7j by using the 
RI-MOM method. It can be also compared with the prediction of one-loop boosted per¬ 
turbation theory Z^®(2 GeV) ~ 0.880. 

We also quote the value of Z^ obtained by using the rough lattice data, without 
implementing the tree-level correction of discretization effects: Z^(2 GeV) = 0.868 ± 
0.003 ± 0.080. The difference in the central value with respect to eq. (ED is less than 
0.5%. As expected, however, the systematic uncertainty is much larger in the latter case, 
due to the signihcantly larger spread of the points in the htting region. In practice, the 
tree-level correction has smoothed the overall behavior of the quark propagator at short 
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distances allowing the rednction of the systematic nncertainty by abont a factor 4, bnt 
affecting the central valne by only a small amonnt. 


5 Extraction of the quark condensate 

One of the advantages of the approach considered in this paper to evalnate the chiral qnark 
condensate is that the renormalization procednre is greatly simplihed: in the OPE of the 
qnark propagator, expressed by eq. o, once the propagator on the l.h.s. is renormalized 
by the qnark held renormalization constant, the r.h.s. tnrns ont to be expressed directly 
in terms of renormalized qnantities. In particnlar, the qnark condensate, renormalized at 
a scale /r, can be extracted directly from the trace of the qnark propagator (i.e. the scalar 
form factor S 2 ) renormalized at the same scale. Fnrthermore, once the qnark propagator 
is improved at 0{a), the operator matrix elements which enter its OPE are antomatically 
improved at the same order. 

In the stndy of the OPE, the physical qnantity which we are interested in is the qnark 
condensate in the chiral limit. To reach this limit, we have followed two procednres. In 
the hrst approach, we extrapolate to the chiral limit the scalar form factor Tj 2 {x) for 
each valne of The qnark condensate is then evalnated by nsing the OPE expressed 
by eq. ©, which is accnrate at the NLO, in the massless case. In this limit, the qnark 
condensate represents the leading term of the expansion. In the second approach, which 
we consider for a consistency check of the caicniation, the order of the extrapolations is 
inverted. At hnite valnes of the qnark mass, the OPE of S 2 at order x^ contains, besides 
the qnark condensate, a term proportional to . In this case, we hrst extract the whole 
0{x‘^) contribntion to the OPE and then extrapolate the resnlt to the chiral limit. As 
we will show in the following, the two procednres yield completely consistent predictions. 
We now discnss the two approaches in more detail. 


Method I: For each valne of the renormalized form factor S 2 (x) is extrapolated to 
the chiral limit, both linearly and qnadratically in either the vector or the axial-vector 
qnark masses. Examples of this chiral extrapolation, for three typical valnes of x"^, are 
shown in hg. H] For each valne of x"^ we have then compnted the qnantity 


Qi{x,n) 


c^^{x,fi) x^/ANc 


+ 0{x‘^) 


and performed a £t to the form 


(19) 


Qi(x,ti) = + Bx^. 


( 20 ) 


Both constant {B = 0) and linear hts have been performed, and the resnlts are presented 
in table 01 see also £g. 01 Since the resnlts of the linear £t are nnstable when the £t is 
limited to the interval {x/aY = [9, 25], we have considered in this case larger distances, 
np to (x/a)^ = 40. In all cases, we find consistent resnlts for the qnark condensate, as can 
be seen from table 01 We also hnd that the contribntion of the 0{x‘^) term is completely 
negligible, and the coefficient B is compatible with zero within the statistical errors. 
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Method I - Chiral Fit 



Figure 4: Quark mass dependence of T, 2 {x) at three different values of x^. The dashed 
lines represent the result of a linear fit. 


Method I - Fit 



Figure 5: Values ofQi{x,p. = 2 GeV) as a function of {x/afi. The solid lines indicate the 
results obtained from a constant fit in x"^ and the estimated systematic error. The dashed 
vertical lines show the range in {x/afi where the constant fit is performed. The result of 
the linear fit in x"^ is also shown (dashed line). 
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= 2 GeV)[ MeV^] 

(x/a)^-range 

Gonstant ht 

Linear ht 


-(265 ±5)^ 

— 

[9,25] 

-(266 ±4)3 

-(265 ±7)3 

[9,40] 


Table 3: Values of the chiral quark condensate in the MS scheme at the scale = 2 GeV 
as obtained from either the constant or the linear fit of eg. 


As a further check of our results, we have also extracted the quark condensate directly 
from the ratio E 2 /S 1 of the two form factors. From eq. (jHl) one hnds that, in the chiral 
limit, this ratio behaves as 




2Nc Cj{x,iJ,) 


+ Oix^). 


( 21 ) 


The determination of the quark condensate from eq. m bypasses the evaluation of the 
quark field renormalization constant Z^. This constant cancels in the ratio, since it enters 
the renormalization of both the form factors Si and S 2 . This also implies that the r.h.s. 
of eq. (ED) is independent of the choice of the renormalization scale fi. We also hnd 
that this ratio, when computed by using the non-corrected form factors, exhibits a more 
stable plateau as a function of The results for the quark condensate obtained with 
the two approaches are in excellent agreement (within less than 2%), indicating that the 
uncertainty connected with the evaluation of is actually negligible. 


Method II: In this second approach we study the OPE of S 2 (a;) at hnite values of the 
quark mass, extract the O(x^) contribution to the expansion and extrapolate it to the 
chiral limit, in order to get the chiral quark condensate. 

The £t of the form factor S 2 to its OPE is shown in fig. |U[ for two values of the quark 
mass. We hnd that the mass term contribution to the OPE, which is leading at very 
short distances where lattice artifacts are more severe, is poorly estimated from the fit. 
For this reason, we have chosen to fix the renormalized quark mass in eq. o to the values 
determined in ref. m and collected in table ^ Therefore, for each value of the quark 
mass we compute the quantity 


Qii{x,m,^) 


S|^^(x, /i) - Cmjx, /i) m^^(p)/47r^ 

C^^{x,^i) x74W 


( 22 ) 


Notice that, in the chiral limit, Qji{x,m, fi) reduces to Qj{x,fi) defined in eq. (IT!I|l . The 
small spread of the results coming from choosing the vector or the axial-vector quark 
masses is included in the systematics. We then £t Qjj{x,m, fi) either to a constant or 
linearly in and extrapolate the result, denoted as p) in £g. El to the chiral limit, 

where it reduces to the chiral quark condensate. The quark mass extrapolation is shown 
in £g. El Though the points in the plot look very well aligned, a quadratic £t in the 
quark mass has been also performed, in order to evaluate the corresponding systematic 
uncertainty. We hnd that the results obtained for the chiral quark condensate with this 
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Method II - X Fit 



Figure 6: Values o/S^(a;,/r = 2 GeV) as a function of (x/a)^ for two different values 
of the quark mass. The solid lines represent the results of the fit to the OPE prediction 
performed by using for the renormalized quark masses the values given in tabled 


Method II - Chiral Fit 



Figure 7: Linear fit of Qifim, fi = 2 GeV) as a function of the quark mass. The extrapo¬ 
lated value is the chiral quark condensate in lattice units. 
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second approach are indistinguishable, within the statistical errors, from those derived by 
using the method I and presented in tab. El 

6 Results and discussion 

Our hnal estimate for the chiral quark condensate is obtained from the results given in 
table El after including the evaluation of the systematic error. We quote: 

= -(265 ± 5 ± 22 MeV)^ , (23) 

where the hrst error is statistical and the second systematic. The latter, which amounts 
to about 25%, is due to: 

- the spread of the points in the fitting regions. As discussed in sect.3, this spread is 
mostly due to discretization effects which are left after the tree-level Cl(a)-correction 
has been applied to the lattice data. This error, of about 18%, represents the main 
source of systematic uncertainty, besides the quenching approximation. 

- Arbitrariness in the choice of the htting interval (within the window a ^ x ^ 
1/Aqcd)- This yields a 4-5% uncertainty. 

- Different functional forms considered in the hts. Performing a quadratic £t instead 
of a linear one in the quark mass extrapolations introduces a systematic difference 
of about 9%. Including the 0{x‘^) contribution in the hts of Qj and Qu to their 
OPEs gives a 4-5% variation in the results. 

- The statistical error associated to the determination of the lattice spacing. This error 
introduces an uncertainty of about 15% in the estimate of the quark condensate. 
Notice that the systematic error associated in the quenched approximation to the 
dependence of the lattice spacing on the physical quantity used to £x the scale is 
not included. We consider this error as a part of the systematic quenching effect. 

- The uncertainty on the quark held renormalization constant Z^, used to renormalize 
the quark propagator. This effect is completely negligible in the determination of 
the quark condensate, as discussed in sect.5. 

- The difference between the results obtained using either the vector or the axial- 
vector definitions of the quark masses. The systematics is slightly affected by this 
effect, by less than 1%. 

- The 10% error on the quenched estimate of Aqcd gives a completely negligible 
uncertainty in the determination of the quark condensate. 

The uncertainty coming from finite volume effects cannot be directly estimated in the 
present study, since our results have been obtained at hxed volume. A study of lattice 
artifacts performed in ref. ra has shown that in the short distance region, which is the one 
relevant for the W-space method, hnite volume effects on the lattice correlation functions 
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in the free theory are negligible with respect to discretization effects. We expect this result 
to remain valid in the interacting theory as well, though a more quantitative conclusion on 
this point would require further investigations. The main source of uncertainty which is 
not evaluated in our estimate of the chiral quark condensate is the effect of the quenching 
approximation. 

In conclusion, in this exploratory study we have investigated on the lattice the OPE 
of the quark propagator at short euclidean distances, and shown the feasibility of this ap¬ 
proach to compute the chiral quark condensate. The result obtained in this way is in good 
agreement with previous determinations of this quantity based on different approaches. 
The strategy investigated in the present study can be also applied to compute on the 
lattice the matrix elements of other local operators which enter the OPE of correlation 
functions at the leading orders. It can be also directly implemented in lattice simulations 
performed with dynamical quarks. 
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Appendix: NLO calculation of the Wilson coefficients 


In this appendix we sketch the NLO QGD calculation of the Wilson coefficients introduced 
in eq. (0. 

The OPE of the quark propagator in euclidean space is expressed by 


T{^p{x)m) 


Cl{x) 7 ^ 7 ^ + —77 Cmix) — 


27r2 

1 

4iV 


C^^{x) + ... 


m\'^ j o 


(24) 


where the dots represent higher powers of and of the quark mass m. All quantities in 
eq. (121 are renormalized at a given scale and in a given renormalization scheme. In the 
following, we will choose the MS renormalization scheme. 

In order to determine the Wilson coefficients at the NLO in QGD, we calculate both the 
left and the right hand side of eq. (12^ up to 0{as) by choosing a common set of external 
states. The coefficients Cj{x) and Cm{x), in particular, can be determined by taking the 
vacuum expectation value of eq. (El in perturbation theory, where the contribution of 
the quark condensate is vanishing. Eq. (El then simply reduces to 


S{x) 


27r2 


Ci{x) 


i^) 


(X 


2\2 


1 ^ 


(25) 
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where S{x) is the quark propagator computed in one-loop perturbation theory. By using 
dimensional regularization, with D = A — 2e, one has 


S{x) = 


d^k 

(27r)^ 


-ik-x 


S(k) 


(26) 


where 


S{k) = Z^- 


L 

I k"^ 


dvr 


7 Zjrn 

Zj.i 


■^'0 






+ 1 


'3-e 


+ 


+ 4 


(27) 


and 1/i = 1/e + log(47r) — From eq. one derives the expressions of the quark 
held and the quark mass renormalization constants in the MS scheme: 


= 1 + ^ Cf\ 
An e 


Z^ = 1 + ^Cf\ 
An e 


By inserting eq. (EH) into eq. EH|l and using 

/■ e-“"’ 1 T{D/2-r) 


(28) 


( 27 r)'^ (k'^Y (An)^/‘^ r(r) 


(29) 


we obtain 


S{x) = 


2n^ 

1 


O^s 

An 

a 


Cni ( 27 i 5 + log(/iV/4)) 




uc 


2\1 


+ 


1 -f Ci. (a - {^-3) (2'Je + 


L dvr 


m 

x"^ 


(30) 


From this result, after comparing with eq. (i2nD, the Wilson coefficients Ci{x) and Cm{x) 
can be readily identihed: 


Qffi 


Ci{x) = Y^fE + log(0xV4) 


a.. 


Cmix) — 1 + —Cp 
47r 


4 - (? - 3) ( 27 e + log(yx^/4)) 


(31) 


In order to compute the Wilson coefficient of the quark condensate, C^^(x), we derive 
a matching equation by inserting both sides of eq. ()24|1 in a connected Green function 
with two external quark helds at hxed momenta: 


((C(^)^J(0))V’7(ffi)V’5(P2)) = 

By putting 


^ {{YiO)YiO))Y^{pi)YsiP2))- (32) 


Cy4^) = Ci4x) + ^C/^^ix) 
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Figure 8: Feynman diagram relevant for the calculation of the Wilson coefficient C^^{x) 
at 0(as)- 


and summing in eq. (El over Dirac {a,ft) and color (a,6) indices, one immediately obtains 
the 0{1) contribution to the Wilson coefficient, 

c?, = 1. (34) 

At 0{as), the matching equation can be schematically written as 


where (.. .)o and (.. .)i represent respectively the 0{1) and 0{as) contributions to the 
Green functions. We now consider the amputated Green functions, and use eq. dSl 
together with the relation 

= I , (36) 

to obtain 

The Feynman diagram which contributes to eq. (nr7|) at 0{as) is shown in hg.|Hl Since the 
matching condition is independent of the choice of the external states, we can evaluate 
this diagram by putting directly pi = P 2 = 0- In addition, by having neglected in 
eq. ipi higher power corrections in the quark mass, we can compute the amputated 
Green functions in eq. (El directly in the limit m = 0. We then hnd: 


= Wtt^CfP 




{27i)D 

{Z^p — 1) — — l) 


^ rf 




fj,i/ 




■I 


+ 

(38) 


where = Z^. By evaluating the Feynman integral with the aid of eq. (j29|l . we hnally 
obtain 

2 - (e + 3) (27 e + log(/i2a;V4))] (39) 


Of C 


Cipp,{x) -1 + ^Cf 


The complete NLO expressions for the Wilson coefficients are derived from eqs. El and 
El by applying the standard NLO evolution functions introduced in eq. dH. 
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